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Abstract. Elemental and isotopic abundances are the fossils of galactic archaeology. 
The observed [X/Fe]-[Fe/H] relations in the Galactic bulge and disk and the mass- 
metallicity relation of galaxies are roughly reproduced with chemodynamical simula- 
tions of galaxies under the standard A-CDM picture and standard stellar physics. The 
isotopic ratios such as ^'^'^^O and ^^'^^Mg may require a refinement of modelling of su- 
pernova and asymptotic giant branch stars. The recent observation of the Carbon-rich 
("1 . damped Lyman a system can be reproduced only with faint core-collapse supernovae. 

Ql4 This suggests that chemical enrichment by the first stars in the first galaxies is driven 

not by pair-instability supernovae but by core-collapse supernovae (~ 20 - SOM©). The 
observed F abundances can be reproduced with the neutrino processes of core-collapse 
supernovae. As in F, the observations of elemental abundances in small systems may 
requires further complications of chemical enrichment. In globular clusters the relative 
contribution from low-mass supernovae is likely to be smaller than in the field, while 
the contribution from massive supernovae seems smaller in dwarf spheroidal galaxies 
than in the solar neighborhood. 

1. Introduction 

Elemental and isotopic abundances are the fossils of galactic archaeology. Different 
elements are produced from stars on different timescales, therefore elemental and iso- 
topic abundance ratios provide independent information on the "age" of a system and 
can be used as a form of "cosmic clock". The formation and evolutionary history of 
galaxies can be constrained in theoretical models by using the observed abundances of 
stars. The space astrometry missions (e.g., GAIA) and large-scale surveys (e.g., HER- 
MES on the AAT) will produce unprecedented information on the ch emodynamical 
structure of the Milky Way Galaxy. Theoretically iKobayashi et all (12006 ) succeeded in 
reproducing the average evolution of major elements (except for Ti) in the solar neigh- 
bourhood including normal Type II Supernovae (SNe II) and hypernovae (HNe). In 
this paper we extend the discussion to isotope ratios (§2) and more realistic chemody- 
namical simulations (§3). From the observed abundances, we discuss the possibility 
of non-standard chemical enrichment in the very metal-poor damped Lyman a system 
(DLA, §4), globular clusters (GCs, §5), and dwarf spheroidal galaxies (dSphs, §6). 

2. Isotope Ratios 

We include the latest che mical enrichment input in to our chemical evolution models. 



which are summarised in IKobaya shi et al.' fe Ollal). The basic equations of galactic 



chemical evolution are described in Kobavash i et all (2006). Figure [T] shows the evo- 
lution of the isotope ratios against [Fe/H] for the solar neighbourhood. In general, 
core-collapse supernovae are the main producers of the major isotopes with more mi- 
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Figure 1 . Evolutio n of isotope ratio s against [Fe/H] for the solar neighbourhood 
with AGB yields. See iKobayashi et all (l2Qlla) for the observational data sources for 
C, Mg, and Ti (dots). The solar ratios are shown with the solar symbols at [Fe/H] 
= in the upper panels. 

nor isotopes synthesized at higher metallicity. For this reason, the evolution of major 
to minor isotope ratios continuously decreases toward higher metallicity. The slope 
changes at [Fe/H] ~ -2.5 and ~ -1.5 are due to the onset of intermediate- and low- 
mass asymptotic giant branch (AGB) stars, respectively. The rapid change in the slope 
at [Fe/H] ~ -1 is caused by Type la Supernovae (SNe la). 

The time evolution of isotope ratios depends on the s t ar forn iation and chemical 
enrichment histories of the system. See iKobayashi et al.l (|2011al ) for the models of 
bulge, thick disk, and halo. In general, the ratios between the major and minor isotopes 
such as ^"^JVig/^^'^^lMg are smaller in the bulge and thick disk, and are larger in the halo 
because of the metallicity effect of supernovae. However, the ^^O/^^O ratio in the halo 
is low due to the production of ^^O from low-mass AGB stars, as also seen in the high 
[(C, F)/Fe] abundances. Therefore, the isotopic ratios can be used as a tool to pick out 
the stars that form in a system with a low chemical enrichment efficiency. This may be 
possible in our halo, but more likely in small satellite galaxies that were accreted onto 
our IVlilky Way Galaxy. 



3. Chemodynamical Simulations 

In a real galaxy, the star formation history is not so simple and the interstellar medium 
(ISM) is not homogeneous. In chemodynamical simulation ( iKobayashi & Nakasatd 
l201lh , hydrodynamics and chemical enrichment are solved self-consistently through 
the galaxy formation. Star formation and chemical enrichment depend on the local 
density and mixing due to dynamical effects such as merging and migration are also 
naturally included. Fig ure [2l shows the frequency distri bution of stars in the solar neigh- 
borhood at present. See IKobayashi & Nakasatd (l201lh for the bulge and the thick disk. 

In our simulations, the metallicity of the first enriched stars reaches [Fe/H] ~ -3. 
At later times, the star forming region becomes denser, and both metal richer and poorer 
stars than [Fe/H] ~ -3 appear. Different from one-zone chemical evolution models, the 
following phenomena occur in the case of inhomogeneous enrichment: i) The age- 
metallicity relation is weak. In other words, the most metal-poor stars are not always 
the oldest stars, ii) SNe la can affect the elemental abundance ratios at [Fe/H] < - 1 
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Figure 2. [X/Fe]-[Fe/H] relations in the solar neighborhood at z = 0. The contours 
show the frequency distrib ution of stars in the sin i ulated galaxies, where red is for 
the highest frequency. See Kobayashi & Nakasatol (l201lh for the observational data 
sources. 

even with the metallicity inhibition of SNe la. The SN la contribution is characterised 
by low [a/Fe] and high [Mn/Fe]. iii) The scatter of elemental abundance ratios becomes 
large if the supernova yield depends on progenitor metallicity such as Na and iv) Some 
of CEMP/NEMP stars can be explained with the local enrichment from AGB stars even 
without the binary effect. In fact, the observed [N/0]-[0/H] trend can be reproduced 
with our simulation without including the effect of rotating massive stars. 



4. Carbon-Enhanced Dampled Lyman a (DLA) System 

The observations of very metal-poor DLAs have opened a new window to study the 
chemical enrichment of the Universe by the first generations of stars. Figure [3] shows 
the elemental abundance ratios from C to Zn relative to Fe. We implement two cases 
with different explosion energies: 1 x 10^^ erg (SN, solid line) and 20 x 10^^ erg 
(HN, sho rt-dashed line) for an initial mass of M = 25Mq and metallicity Z = (see 



Kobayashi et al. 2011b for the details). An efficient mixing-fallback is adopted in both 



models, black hole masses are as large a s - 6Mc:^, and thus the ejected iron mass is 
much smaller than "normal" SNe/HNe in 'Kobaya shi et al l ('2006') that are responsible 
for the chemical enrichment in the Galaxy. Since C is synthesized in the outermost 
region of the ejecta, [C/Fe] ratio is higher for faint SNe/HNe than "normal" SNe/HNe. 

We then compare these model calculations wit h the availabl e abundance measure- 
ments in the extremely metal-poor DLA reported by lCooke et al.l (HoiOb, filled circles). 
This DLA was originally identified in the SDSS spectrum of the QSO J0035-0918. 
Follow-up high resolution spectroscopy indicated Zabs = 2.3400972, log A/^(HI) /cm"^ = 
20.55 + 0.1, and [Fe/H] ^ -3.04. 1) A pronounced carbon enhancement [C/Fe] ^ +1.53 
is found, which is well reproduced with our faint SN/HN models. Low-mass AGB stars 
(1 - 4Mo) can also provide such high [C/Fe]. However, such low-mass stars are un- 
likely to contribute at the redshift of the C-rich DLA. With the star formation history in 
the solar nei ghborhood, [C/Fej reach es the maximum value at z = 1.8 due to the AGB 
contribution ('Ko bayashi et aklEoilah . 2) The [O/Fe] ratio is higher than for non-faint 
supernovae ([O/Fe] ~ 0.5 - 0.6 for 25Mo) and is consistent with the faint SN model 
because of the smaller ejected Fe mass. 3) The [(Si,S)/Fe] ratios are similar to those 
of non-faint supernovae, and also consistent with the faint SN/HN models. 4) The low 
Al abundance strongly suggests that the enrichment source is not Pop II supernovae 
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Figure 3. The elemental abundance pattern of the metal-poor C-rich DLA (filled 
circles) and peculiar DLA (open circles). The solid and short-dashed lines show 
the nucleosynthesis yields of faint core-collapse supernovae from 25 M© stars with 
mixing-fallback. The dotted line is for PISNe from ITOM© stars. 

but primordial supernovae. 5) The low N abundance is also consistent with the faint 
SN/HN models, and the high [C/N] ratio cannot be explained with mass loss from ro- 
tating massive stars or intermediate-mass AGB stars (>4Mq). 

In order to discuss the detailed explosion mechanism of the supernova, it is nec- 
essary to obtain the elemental abundances of iron-peak elements. For the C-rich DLA, 
because of the low metallicity, it is impossible to detect heavier elements t han S except 
for Fe. In Fig. [3] we overplot the element abundances recently reported by lCooke et al. 



(I2010al . open circles) for the Zabs = 1.62650 DLA in front of the gravitationally lensed 
quasar UM637A. This abundance pattern can also be explained with faint supernovae. 

The dotted line is fo r the nucleosynthesis yie lds of the pair instability supernova 
(PISN) of a 170Mo star (lUmeda & Nomotol[2002h . Compared with core-collapse su- 
pernovae, the abundance pattern of PISNe can be summarized as follow: 1) The odd-Z 
effect is much larger than ~ 1 dex. 2) [(Si, S, Ar, Ca)/Fe] are much larger than [(O, 
Mg)/Fe] because of more extensive explosive oxygen burning. 3) [Cr/Fe] is much 
larger because of the larger incomplete Si-burning region. 4) [(Co, Zn)/Fe] are much 
smaller because of the much larger ratio between the complete and incomplete Si- 
burning regions. All of these characteristics disagree with the observed elemental abun- 
dances of these metal-poor DLAs. Even at high-redshift, there is no signature of the 
existence of PISNe. The [Si/C] for PISNe is as large as +1.5, which is also inconsis- 
tent with the observational estimate in the intergalactic medium (IGM, [Si/C] ~ 0.77). 
The IGM abundance looks more consistent with normal (non-faint) core-collapse su- 
pernovae with [C/Fe] ~ and [Si/Fe] ~ 0.7 (Kobayashi et al. 2006). 

In the early stages of chemical enrichment, the ISM is supposed to be highly inho- 
mogeneous, so that the properties of the first objects can be directly extracted from the 
comparison between the observed elemental abundances and nucleosynthesis yields. 
Since the DLA abundances reflect the chemical enrichment in gas-phase, the binary or 
accretion scenarios of the BMP stars do not work. To explain the abundance pattern of 
the C-rich DLA, enrichment by primordial supernovae is the best solution. It is inter- 
esting that the observed DLA abundance is very similar to those of BMP stars in the 
solar neighborhood including the ultra metal-poor star HB0557-4840 ([Fe/H] = -4.75, 
[C/Fe] = +1.6) and BD+44°493. Some of BMP stars in dSphs and the Galactic outer 
halo also show similar carbon enhancement at [Fe/H]^ - 3. Chemical enrichment by 
the first stars in the first galaxies is likely to be driven by core-collapse supernovae from 
~ 20 - 50Mo stars, which may be co nsistent with the latest simulations of primordial 
star formation (e.g. JGreif et aLl[20TTh . 
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Figure 4. Evolution of the [F/0] ratio against [0/H] for the solar neighbourhood 
with SNe II, HNe, and SNe la only (short-dashed lines), with AGB stars (long- 
dashed lines), with the v-process of SNe II and HNe (solid line and dot-dashed line 
for = 3 X 10^^ and 9x10^^ er g, respectively). The dotted line is for the model for 
GCs. See lKobayashi et al] (|2011q) for the observational data sources. 

5. Fluorine Problem and Globular Clusters (GCs) 

Fluorine is an intriguing, though currently poorly studied element. ^^F is mainly pro- 
duced by core and shell He-burning at T > 1.5 x 10^ K in low-mass and massive stars. 
The production of F in AGB stars is highly mass dependent, where F production peaks 
at ~ 3M0 at solar metallicity; in higher mass models F is destroye d by a-captures 
caused by the higher temperatures reached during He-burning. In i Kobayashi et all 
(l2011ch . we have calculated the nucleosynthesis yields of core-collapse supemovae 
(SNe and HNe) including the y-process. IMost of the energy from core-collapse super- 
novae is released as neutrinos and anti-neutrinos (> 10^^ erg). The y-process does not 
affect the yields of major elements such as Fe and a elements, but it increases those of 
some elements such as B, F, K, Sc, V, IMn, and Ti. 

Figure |4] shows the evolution of [F/O] against [O/H]. Without the AGB yields 
and the y-process (short-dashed line), the predicted F abundance is too low to meet 
the observational data at all metallicities. With the AGB yields (long-dashed line), 
[F/O] shows a rapid increase from [O/H] > - 1.2 toward higher metallicities, which 
corresponds to the timescale of 2 - AMq stars in the solar neighborhood. At [O/H] ~ 0, 
[F/O] reaches -0.14, which is 0.26 dex larger than the case without the AGB yields. 
However, the present [F/O] ratio is still significantly lower than the observations at 
[O/H] ~ 0. Note that the F yields from AGB stars were increased with updated reaction 
rates. 

The timescale of supernovae is much shorter than AGB stars, which means that 
the [F/O] ratio at low metallicities can be strongly enhanced by the y-process occurring 
in core-collapse supernovae. With the standard neutrino luminosity £'y = 3 x 10^^ erg 

(solid line), the [F/O] ratio shows a plateau of [F/O] 0.4 at [O/H] < - 1.2, and 

reaches [F/O] ~ +0.19 at [O/H]^0. This is consistent with the observational data of 
field stars at -0.5 < [O/H] < 0. If we adopt a larger Ey - 9 x 10^^ erg (dot-dashed 
line), [F/O] can be as large as ~ +0.37 at [O/H] ~ 0. In the bulge the star formation 
timescale is shorter and the average metallicity is higher than the solar neighborhood, 
but the [F/O] ratio is not so different at [O/H] ~ 0. The observations for the bulge stars 
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(filled circles) might suggest that the initial mass function is also different, although the 
number of observations is too small to make a conclusion. 

At -1^ [F/0]< - 0.5 the observational data are for stars in GCs, which seem to 
be more consistent with the models with the AGB yields only than with the y-process. 
However, it is unlikely that the existence of the y-process depends on the environment or 
metallicity. The neutrino luminosity may be small in the case of faint supernovae with 
a large black hole, which give high [a/Fe], but there is no significant difference seen 
in the [a/Fe] ratio between field halo stars and GC stars. One possible scenario is as 
follows: in GCs, the contribution from low-mass supernovae is smaller than in the filed. 
Since the star formation occurs in a baryon dominated cloud with very high density, the 
initial star burst can be very intense. After the initial star burst, because of the small 
gravitational potential, galactic winds are generated immediately after the explosion of 
massive supernovae, which may remove the contribution from low-mass supernovae. 
The small production of a elements from low mass supernovae means that the [F/O] 
ratio can reach values as large as ~ 0. In contrast, massive supernovae produce more a 
elements which results in [F/O] ratios of ~ -0.5, consistent with the observational d ata. 
The dotted line shows an example of such a GC model (see lKobayashi et aDl201ld for 
the details). 

At low metallicities ([O/H]^ - 1.2) F is produced only from supernovae, and thus 
the observations of field stars at low-metallicities are important for constraining the 
neutrino luminosity released from a core-collapse supernova. If the neutrino luminosity 
is specified, the F abundance along with C could be a good clock in the study of galactic 
archaeology to distinguish the contribution from AGB stars and supernovae. 

6. Dwarf Spheroidal Galaxies (dSphs) 

The F observations of stars in GCs suggest that the star formation and chemical enrich- 
ment histories of GCs are different from those of field stars and that low-mass super- 
nova played a smaller role in shaping the chemical evolution of these systems. This 
is the opposite to th e situation for dSphs. DSphs show low [a/Fe] and low [Mn/Fe] 
(iRomano et al.ll201ll ). of which elemental abundance pattern is more c onsistent with the 



enrichment from low-mass SNe II than that from the SN la enrichment (|Kobayashi et al.l 



l2006l . see also Figure 20 of lKobayashi & Nakasatoll20TTl ). There is no strong metallicity 



dependence on Mn yields of SNe la (T. Ohlubo, K. Nomoto et al, in prep.). In dSphs, 
the dark matter component is large, the gas density is low, the star formation rate is low, 
and thus the contribution from massive supernovae is expected to be smaller than in the 
Milky Way halo. The ele mental abundance p attern (low [(Cu, Zn)/Fe]) of anomaly stars 



in the Galactic halo (Nissen & Shusteijl2011h is also consistent with low-mass SNe II. 
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